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Abstract: We investigate the breakdown of collinear factorization for non-inclusive ob- 
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Coulomb/Glauber gluon exchanges with other soft exchanges. Beyond QCD, factorization 
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1 Introduction 

Factorization theorems are a cornerstone of QCD calculations. In particular, it is the 
inclusive collinear factorization theorem that allows the calculation of cross sections in 
hadron-hadron collisions as the convolution of universal (i.e. process independent) parton 
distribution functions and (perturbatively calculable) partonic cross sections. 

The factorization theorem has been proven for sufficiently inclusive observables in the 
final state of the scattering of colourless hadrons [1]. However, it is often assumed that 
partonic scattering amplitudes also factorize, e.g. in the calculation of amplitudes to high 
orders and their resummation to all orders. Moreover, this assumed factorization forms the 
basis upon which parton shower Monte Carlo event generators are constructed. In neither 
of these cases is the factorization guaranteed. Indeed, in a recent paper, Catani, de Flo- 
rian and Rodrigo (CdFR) explicitly wrote down amplitude and cross section contributions 
that violate collinear factorization for non-inclusive processes with two incoming coloured 
partons [2]. 

In this paper we aim to shed light on the breakdown of collinear factorization. To 
that end, we re-consider the results presented in CdFR. Most of their paper concerns the 
question of whether there are factorization-violating terms in higher-order corrections to 
QCD amplitudes. The answer, for initial-state (space-like, SL) collinear branching in a 
process with two incoming coloured partons, is definitely Yes. CdFR derived an all-orders 
expression for the singularities of the collinear splitting amplitude and gave its explicit 
realization at one- and two- loop order. We discuss the structure of these results and isolate 
the factorization violating terms in a way which makes it clear that their physical origin is 
the Coulomb interaction between the two incoming partons long before the hard interaction 
or between two outgoing partons much later than the hard interaction. 

CdFR then discuss whether these factorization violations at amplitude level lead to 
factorization violation at cross section level. The answer at one-loop order is definitely No. 
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CdFR does however calculate a contribution that violates cross section factorization at two- 
loop order. We show below that this result is, in fact, only non-zero for certain electroweak 
processes. For QCD scattering processes, factorization violation effects cancel and the 
answer at two-loop order remains No. In the following, we argue that this will not continue 
to be the case at still higher orders, and we do expect factorization violation at three-loop 
order in the SL collinear limit in parton-parton cross sections. We specify the form we 
expect for these terms. 

This paper is set out as follows. In the remainder of this introduction, we briefly 
recall our work on the logarithmic structure of high-order corrections to non-inclusive QCD 
observables such as dijet production with a jet veto (often referred to as 'gaps between jets') 
[3, 4]. As we shall see, this is relevant because the "super-leading logarithms" we identified 
there have the identical structure and physical origin as the factorization violating terms 
discussed here. Then, in section 2, we recap some of the formalism used by CdFR in 
deriving results in the SL collinear limit. In section 3 we review the results of CdFR for 
the one-loop corrections to factorization of SL collinear emission at amplitude level. We 
rewrite the final result in a form in which it is manifest that the physical origin of the 
factorization violating term arises from an effective Coulomb interaction between the two 
incoming partons. We show that it does not contribute to physical cross sections at this 
order. In section 4 we consider the two-loop corrections at amplitude level. We confirm 
the CdFR results and again rewrite them in a form in which the connection with Coulomb 
gluons is manifest. In section 5 we show that the factorization violating terms do not 
contribute to QCD cross sections at this order either. We discuss the conditions under 
which they could give a non-zero contribution and, within the Standard Model, show that 
these are only satisfied by certain electroweak processes. Finally, in section 6 we make 
some concluding remarks, including a prediction of the form of the term that will violate 
the cross section factorization of collinear singularities at three-loop order (for coloured 
incoming partons). 

In Refs. [3, 4] we studied the structure of the 'gaps-between-jets' cross section in QCD. 
As is well known (see for example the series of papers by Sterman and collaborators [5-7]), 
"leading" logarithmic terms ~ a™ log n (Q/Qo)) where Q is the hard scale and Qo the veto 
scale, can be summed to all orders by a matrix-valued evolution in colour space. The 
anomalous dimension matrix that drives this receives a real (strictly, Hermitian) contri- 
bution from virtual gluons with light-like momenta and in what follows we shall refer to 
these as 'eikonal gluons'. It also receives an imaginary (strictly, anti-Hermitian) contri- 
bution from space-like virtual gluons that correspond to rescattering between the partons 
long before or after the hard scattering process. We refer to these as 'Coulomb gluons', 
although they are also referred to as Glauber gluons in the literature. 

In Ref. [3], we considered the first correction to this picture. Specifically, we allowed for 
one additional gluon, either real or virtual, with momentum corresponding to phase space 
regions lying outside the jet veto region 1 . Following the work of Dasgupta and Salam, [8], 

lr rhe virtual corrections in the anomalous dimension matrix are instead integrated over the region of 
phase space where real emissions are forbidden, i.e. inside the jet veto region. 
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we expected to find a tower of additional leading logarithms, called non-global logarithms. 
Physically these arise from the requirement that although radiation outside the veto region 
is allowed, secondary radiation from it back into the veto region is not. We did find these 
non-global logarithms. To our surprise we also found that, even when the 'out-of-gap' gluon 
is far forward (i.e. collinear to an incoming parton), the subsequent colour evolution differs, 
depending on whether the low-angle gluon is real or virtual. This triggers a breakdown of 
the 'plus prescription' and to the conclusion that collinear logarithms can only be factorized 
into the incoming parton density functions at or below the scale Qq. In subsequent confer- 
ence presentations [9], we described this as a breakdown of QCD coherence and discussed 
further the fact that it signals a breakdown of collinear factorization for processes with in- 
coming coloured partons. In the gaps-between-jets calculation, this mis-cancellation means 
that the integration over the momentum of the out-of-gap gluon is double-logarithmic and 
leads to the appearance of super-leading logarithms: one out-of-gap gluon contributes a 
tower ~ a™log n+1 (Q/Qo), n > 4, and we speculated that the leading behaviour at high 
order is actually ~ a™ 

log 2n " 3 (Q/Q )- This speculation was subsequently confirmed, at 
n = 5 with two out-of-gap gluons, in a fixed-order calculation [10]. Ref. [4] re-performed 
the calculation using the colour basis independent notation [11] also used by CdFR and 
showed that it is a general feature of QCD scattering amplitudes. These effects were also 
considered for event shape variables in Ref. [12] where they were called "coherence violat- 
ing logarithms" since, depending on the class of event shapes, they can be super-leading, 
leading or sub-leading. Importantly, these effects are generic to all non-inclusive final state 
observables (and they are certainly not unique to non-global observables) . 

As emphasized in Ref. [4], the origin of the mis-cancellation between the evolution of 
a system with a real or virtual collinear gluon comes entirely from a mismatch between the 
colour matrices of the Coulomb gluon contributions to the two amplitudes. For inclusive 
observables, this would only lead to a phase (a purely anti-Hermitian transformation in 
colour space) in the amplitude and would not contribute to physical cross sections. How- 
ever, the non-Abelian nature of QCD means that the colour matrix of the Coulomb gluon 
contribution does not commute with that of other soft (eikonal) gluon exchanges that gen- 
erate the Hermitian part of the evolution in the case of non-inclusive observables. This 
non-commutation can (and does) lead to Coulomb-induced factorization violation in phys- 
ical cross sections. This occurs first at fourth order relative to the hard process, since we 
need one collinear emission, two Coulomb exchanges (to give a real Hermitian contribution, 
~ (i7r) 2 ) and one real exchange (for them to not commute with). 

In the following, we will show that the physical origin of the non-factorizing terms 
identified by CdFR is identical to that of the super-leading logarithms and that the same 
argument about their leading effect on physical cross sections applies: they can contribute 
to fourth-order corrections to hard QCD processes, i.e. three-loop corrections to the single- 
collinear splitting functions. 
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2 Space-like collinear factorization 



At tree level, the factorization theorem for QCD amplitudes in the multiparton collinear 
limit reads 

A*(°>\«Sp(°> M {0) ). (2.1) 



is the tree-level amplitude 2 for some process involving n coloured partons. We follow 
the convention of CdFR in defining all momenta as outgoing, so incoming partons have 
p® < 0. To be concrete, when discussing processes with two incoming partons (as we are 
throughout this paper) we assume that partons 1 and m+1 are incoming and all others are 
outgoing. [A^ )) is the representation of Ai^ in colour space [11]. Equation (2.1) defines 
the factorization theorem of the amplitude Ai^ in the multicohinear limit in which m 
of the partons become collinear (in the SL case of interest, one of the m is the incoming 
parton 1). The accuracy is as specified in CdFR: in the multicohinear case it gives the 
dominant contribution in which all scalar products among the m collinear partons are 
of the same order and all vanish together, and does not include sub-dominant (but still 
singular) contributions. In the particular case of the two-parton collinear limit there are no 
sub-dominant corrections and this formula captures the entire singular collinear behaviour. 
A4^ is the tree-level amplitude for a process in which the n — m non-collinear partons' 
momenta and colours are the same as in Ai^°\ but the m collinear partons are replaced 
by a single (on-shell) parton with momentum 

V 2 m 



fjet 

,2 



t=l 



where n is a light-like (n = 0) vector introduced to define the collinear limit. In the SL 
case of interest, P° < 0. Flavour conservation in the collinear limit defines the flavour of 
parton P and colour conservation implies that it has colour 

m n 
i=l j=m+l 

Note that in general, colour conservation means that, when acting on physical states, the 
colour operators obey 



]TT 4 = 0, (2.4) 



i=l 

from which the second equality of Eq. (2.3) follows. 

Sp^ is an operator describing the collinear splitting P(Tp) — > pi(T\) +p2(T2) + ■ ■ • + 
p m (T m ). It can be represented by a non-square matrix in colour space, since it acts on the 
space of n — m+1 partons to produce a state in the space of n partons. Equation (2.1) is 
defined as "strict" factorization, because Sp^ depends only on the colours and momenta 



2 Strictly speaking, we mean the lowest order non-zero amplitude, which for some processes could mean 
a loop amplitude, for example gg — > "fg. 



- 4 - 



of the m collinear partons and depends only on the colours and momenta of the 

non-collinear partons and P. 

CdFR show that the structure of Eq. (2.1) continues to all orders: 



M) ~ S P 



M), 



(2.5) 



where the kets \M) and \M) and the operator Sp have perturbative (loop) expansions 
starting from loops, i.e. 



M^) = Y J Sp {V 



l'=0 



(2.6) 



Equation (2.5) is a factorization, but it is described as a "generalized" factorization, because 
although Ai still depends only on the colours and momenta of the non-collinear partons 
and P, Sp depends in general on both the collinear and non-collinear partons. 

In this paper we study the infrared poles (i.e. e poles in dimensional regularization) 
of the multi-collinear splitting matrix Sp at various orders 3 . In fact, the only ingredients 
that we need to do this are Eq. (2.5) and the results for the singularities of the on-shell 
amplitudes |At^) and \M )■ These have been known at 1 loop (I = 1) for many years, 
and were first written down in the colour basis independent notation in Ref. [11]. They 
were written down for the two-loop case in Ref. [13] and have recently been explored at 
higher loops [14-18]. We define the singular parts of the Moop amplitude as 



M il) ) = J2 lin M {l ~ V) )+ M (/)fin ' 



(2.7) 



with analogous definitions for operators 1 1 ^ for the reduced amplitude We give 

explicit expressions for these operators in the appropriate sections below. 

We can use the factorization of amplitudes in the collinear limit to derive the factor- 
ization of cross sections: 



M 



M 



M 



P = SpT Sp 



(2.8) 



If Sp obeys strict factorization, then P does too. On the contrary, it is possible that strict 
factorization is broken at the amplitude level, but nevertheless the factorization-violating 
terms do not contribute to physical cross sections. In fact this is true for time-like (TL) 
collinear splitting, less trivially for SL collinear splitting in processes with only one incoming 
coloured parton (e.g. deep inelastic scattering) and, less trivially still, for general processes 
at one-loop level and QCD processes at two-loop level, as we will show. 

We close this section by noting a property of QCD tree- level matrix elements \A4 ° ) 
(and also 5p (0) ) that we will make use of in the following. In the colour basis independent 



3 Of course, taking m = 2 we recover the two-parton collinear limit. As shown in CdFR, in this case it 
is possible to calculate the non-factorizing terms exactly in d dimensions but, while important for detailed 
calculations, we do not believe that this illustrates any additional physics, so we work in the multicollinear 
case throughout. 
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notation, one works in an orthonormal basis for the colour states of a given scattering 
process. It is natural (and has been done in all calculations to date) to define these 
basis states as combinations of generators, delta functions and structure constants (if abc 
and d abc ) with real coefficients. It was noticed in Ref. [19] that, in these natural bases, 
all anomalous dimension matrices calculated in the literature are (complex) symmetric 
matrices. In Ref. [20] it was proved that this is because the matrix representation of any 
colour operator Ti-Tj (which is Hermitian in any orthonormal basis) is real in these bases, 
and hence symmetric. In the following we will exploit the existence of such a basis, even if 
it is not necessary that it be specified explicitly. 

The vector representation of the amplitude |.M^) is given by its projections onto 
the basis vectors of the space. An important point in the following is the fact that tree- 
level QCD amplitudes do not introduce any additional phases into this representation. 
That is, all the elements of the vector representation of have the same phase and 

so the outer product | A4 ) { A4 ^ | (which is called the "hard matrix" in the language 
of Sterman et al.) can be represented as a matrix with real entries. Since it is also 
Hermitian it must be symmetric. The same argument applies to Sp^°\ It is a tree-level 
QCD operator whose matrix representation is a transformation between two orthonormal 
spaces (of different dimensions). If both of those spaces are represented by natural bases, 
the matrix representation of Sp^ is real. 

3 Space-like collinear factorization at one loop 

At one loop [11], the operator i" can be written 4 



2vr 2 ^ V e 2 * e "I e ^ ' J V a 



2 



(3.1) 



with an exactly analogous expression for J 1 . Examining the analytic structure of the 
logarithm, we see that when i and j are both incoming partons, or both outgoing partons, it 
contributes an imaginary part. Physically, this imaginary part corresponds to an absorptive 
contribution coming from on-shell parton-parton scattering long before, or long after, the 
hard scattering process. One can rewrite this equation to make the imaginary part manifest, 
for example as 



2ir 2 

where 



t (? c < + h + t c ') - 7 1 T - • T = h (^) + 2J f T - • («) 



Ti+T m+1 = - T i ( 3 - 3 ) 

i^l,m+l 



4 To avoid a cluttered notation, we do not explicitly indicate that J, and all the related operators we 
consider, are renormalized and only have poles of IR origin. Thus we drop the R superscript used by CdFR, 
as well as the argument of as, which is the MS coupling evaluated at /i . Finally, we also drop the argument 
of /(e), except where it is not e, when we explicitly write it. 
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is the total colour in the s channel. In this form, one can see that the contribution comes 
symmetrically from initial-state and final-state interactions, which, since colour is con- 
served, are equal. It is more useful for practical calculation to write this entirely in terms 
of the initial-state contribution. Moreover, purely Abelian imaginary parts ~ mCi can- 
not contribute to cross sections. Therefore, in the following, we shall write the physically 
important imaginary part as ^-T\ ■ T m+ i. 

Writing the operator as in Eq. (3.2), we are separating different regions of the inte- 
gration over gluon loop momenta: "eikonal" gluons (i.e. with on-shell momenta) give rise 
to the logarithm, whilst "Coulomb" gluons (i.e. with space-like momenta) correspond to 
on-shell external parton-parton scattering and give rise to the in term. It will turn out 
that the factorization breaking contributions can be understood entirely in terms of the 
non-trivial colour operators associated with Coulomb and eikonal gluon exchanges. 

The singularities of the one-loop splitting matrix, defined by 



Sp 



« = lg> Sp(°> + Sp^ 



(3.4) 



can now be extracted using 



iffs P ^ = iWsp^-s P m {1) 



(3.5) 



which we write as 5 



r (i) 

c 



j(D 



r(!) 



(3.6) 



After some colour algebra, including the colour conservation condition in Eq. (2.3), we 
obtain 



c 2vr 2 



7i c p + hp 




i=2 



with 



•(1) 



1 m 

- Y Ti-T e \n 



\Si£\ 



+ A 



(i) 
c ' 



a 



■ 2 x — Twh • (Ti 

Z7T e 



(3.7) 



(3.8) 



All notation not explicitly defined here is taken over directly from CdFR, but in brief: C, 
is the Casimir of particle type i, 7$ is a flavour dependent real constant and Zi = pi-n/P-n 
is positive for the incoming collinear parton (i = 1) and negative for the outgoing collinear 
partons {i = 2 . . . m) . 

Note that Eq. (3.7) is exactly equal to CdFR Eq. (5.31), but it is rewritten in such a 
way that our is different to their A^, in Eq. (5.32). The terms explicitly written in 
Eq. (3.7) obey strict factorization: they depend only on the kinematics and colours of the 
collinear partons. On the other hand, A^ violates strict factorization since it depends on 



3 See the discussion above Eq. (5.49) in CdFR. 
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the colour of a non-collinear parton, i.e. the other incoming parton, m + 1. Written in this 
form, it is very clear that the violation of strict factorization observed by CdFR is directly 
related to a mismatch between the colour structures of the Coulomb gluon contributions 
in the full (T^i • T±) and factorized (T^i • Tp) matrix elements. All other terms obey 
strict factorization. This mismatch was illustrated diagrammatically in Figure 1 of Ref. [4] 
and in the associated discussion. 

Turning to the factorization theorem for cross sections, we have 

P« = flfp(°)t Sp« + Sp^ Sp® = Sp™ (lg> + /g )f ) SpM + pd)fi- (3.9) 
= Sp(°)t/«Sp(o)+p(i)^, (3.10) 

1 1 \ m ( 1 1 2 \ 1 m 

72 C P + - 7j5 ) - £ Ud + - 7 , - -Qlnh| --E T i' T < 



p 2vr 



In- 



(3.11) 

We see that the divergent terms do now have a factorized structure. This is because 

the non-factorizing term in Iq , is anti-Hermitian and it has cancelled in the combination 

Iq + Iq ■ The only possible source of factorization violation is in the finite term, 

P (i)fin. = Sp m S pW &n - + h.c , (3.12) 

since in the multi-parton collinear limit it is not (yet) proven that Sp^ &n ' contains only 
(finite) factorizing terms. In the two-parton limit on the other hand, Sp^ &n ' has the same 
structure as Sp^ and hence factorizes. 

In summary, in this section, we have presented results that completely agree with 
those of CdFR, but we believe that our writing of the non-factorizing contribution to 1^ , 
Eq. (3.8), displays its physical origin more clearly. 



4 Space-like collinear factorization at two loops: amplitude level 

At two loops [13], the operator i" is 



j(2) 



J(D 



a s 1, 
+ 2i +7"° 



}9 1 n 
+(£) 7£^ 2) < c 4 - 1 ) 
i=l 



(2) (2) 

with an exactly analogous expression for I . K and H) are more (real) constants that 
may be found in CdFR. Recall that we suppress the e argument of J®, except where its 
argument is not e. 

The singularities of the two-loop splitting operator can then be extracted as 



Sp 



(2) 



+ finite terms, 



where 



^ (2)di ' 



j(l) )5p (l)fin. 



(4.2) 
(4.3) 
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which vanishes if there are no factorization violating terms in Sp 



(l)fin. 



(2) 

We can extract I c as 



r(2) 

-c 



i.e. 

r(2) 
1 C 



1 

~2 
1 

+ 2 



r (D 

■c 



r(l) r(l) 
1 > J C 




/(2) 



jg)( 2e )-J^ +KI${2e) + 



7 (2) +7 (1) (/« 



r(l) 



(4.4) 



'2,1 + «• I 
2vr e 



E^ 2) -4 2) 



(4.5) 



Note that all of the parts that are non-trivial colour operators are determined by one-loop 
operators. Thus, at two-loop order, the violations of factorization can be understood in 
terms of one-loop violations. To this end, it is useful to write 



r (i) 

■o 



r (l)fact. 

C 



+ A 



(1) 

c ■ 



(4.6) 



r (l)fact. 

C 



depends only on the collinear partons' colours and momenta and has both Hermitian 



and anti-Hermitian parts, while A^ depends on the colours of collinear and non-collinear 
partons and is purely anti-Hermitian. We then have 



r (2)fact. 




jg)fact. (2e) 



r (l)fact. 

-c 



+ KI 



(l)fact. 
C 



(2e) 



(2) 

and a factorization violating contribution to I c of 



L (2) 



r (2) 
c 

1 

~2 



r (2)fect. 

-c 

r (l)fact.r(l) 



■(1) r (l)fact. 



+ 2 



1 Tt(i: 



A {1) 



a s ( 1 , 



Ag } (26) 



(4.7) 



(4- 



(4.9) 



Note that, unlike A^ , which is purely anti-Hermitian, A^ contains both Hermitian and 
anti-Hermitian parts. 

In summary, in this section, we have shown results that completely agree with those of 
CdFR, but again they are slightly rewritten to isolate more clearly the physical origin of the 

(2) 

non- factorizing contributions to I c : they arise entirely from non- factorizing contributions 
that are already present in the one-loop results. 



5 Space-like collinear factorization at two loops: cross section level 

In order to compute cross sections, we need 

P< 2 > = 5pW Sp&> + Sp«t Sp« + SpW 5p(°) 



,(2) 



(2) 

+ f + finite terms, 



(5.1) 
(5.2) 
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with a factorizing contribution of 



pf = s P w /g )fact - + 



r (l)fact. 
1 C 



r (2)fact.t 



r (l)fact.f 
1 C 



2 



+ jg)fect.f jg)£act. I g p (p) 

(5.3) 



Any factorization violation at two-loop order arises from the remainder: 

,(2) = o„(o)t / a I 2 ) + A(?)t , rXWV _ aW r« fact -t j_ rW fact -AW L ««(°) 



+5p (0)t {/^5p (1)fin - + 5^ (2)div -} + h.c. + finite terms. (5.4) 

The first line contains all of the IR divergences associated with divergent operators acting 
on the tree-level splitting operator, while the second line contains possibly non-factorizing 
divergences due to non-factorizing finite parts of the one-loop splitting operator (recall that 
5^ (2)div ' depends on Sp^ &n -). 

Inserting the expression in Eq. (4.9) into the first line of Eq. (5.4), we obtain 



P ( 2 ) = Is p (0)t + jWT + jrgJM. + jg^ct-t j > £g>j Sp<°> + second line. (5.5) 

This important expression has a very simple physical interpretation. It expresses the fact 
that Coulomb exchange (which is the physics of A^) does not commute with eikonal 
exchanges between any two partons in the matrix element (which is the only relevant 
physics of 1^ and since the Coulomb exchanges in these operators cancel after adding 
the Hermitian conjugate 6 ). 

We can relate this commutator to that which is explicit in the derivation of super- 
leading logarithms in Ref. [4] by considering the limit in which the partons of the reduced 
matrix element divide into two clusters that are separated by a large interval in rapidity, Y. 
In this case, the only large s« are associated with the case where partons i and j are on 
opposite sides of the interval. Taking the leading Y behaviour is therefore equivalent to 
neglecting 1^ and simplifying 

T (i) +J (i)t^2 y2l) 
2tt e 

where T% is the colour exchanged across the rapidity interval. Equation (5.5) then simplifies 
to 

Pg oc iirYSpW [T?, TVh • (T l - Tp)] Sp® . (5.7) 

This equation embodies exactly the same physical result as Eq. (4.13) of Ref. [4]: the 
violation of factorization is driven by the non-vanishing commutator of the eikonal gluon 
contribution with the mismatch between the Coulomb gluon contributions of the full and 
factorized matrix elements. To make this more precise, we should explain the connection 
between the notation used here and that used in Ref. [4]. In that paper, the operator 
is the colour part of Sp^°\ Ti • T2 is • T\ and ti • t2 is T^x ■ Tp. Finally, since 



r(l) . rWt 1 r W fact - , r(l)fact.A A (l) 



'For that reason the superscript 'fact.' is actually redundant in Eq. (5.5). 
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Tm\-\ ■ Tp and T\ commute with Sp^ (in the sense defined more precisely in Eq. (4.12) 
of Ref. [4]), Eq. (5.7) can be written identically to Eq. (4.13) of Ref. [4]. 

In general, the commutator in Eq. (5.5) is non-zero. However, we can exploit its 
properties as a matrix in colour space to show that its expectation value on a QCD tree 
amplitude Sp (0) |7W ( ') is zero. As discussed at the end of section 2, there exists a basis 
in which the matrix representation of any colour operator Tj • Tj is real. Therefore, in 
such a basis, the commutator in Eq. (5.5) is purely imaginary. Since it is also Hermitian, 
it must be antisymmetric. As we have noted, in such a basis QCD tree amplitudes are 
real and hence |jM^) (A4^°^| is symmetric. The final step of the argument is to note that 
(Ai^\A\A4^) can be written as Tr (|7W (0) ) (A4 (0) |)A and 



Tr 



SA 



(5.8) 



for any symmetric, S, and antisymmetric, A, matrices. Since there exists a basis in which 
the result is zero it must be zero in any basis. 

Therefore, although the commutator in Eq. (5.5) is non-zero, its expectation value on 
a QCD tree amplitude is zero and 

{ iy\ 

P^f = second line (pure QCD processes). (5-9) 

We recall that the terms in the "second line" of Eq. (5.4) arise from divergent one-loop cor- 
rections to finite one-loop corrections to the splitting operator, which are not yet proven 
to factorize. We have shown that, apart from this possible source, factorization violat- 
ing terms in the one- and two-loop splitting amplitudes do not contribute to factorization 
violation of QCD cross sections at the one- or two-loop level. 

CdFR considered the most general case, in which it is true that the expectation value 
of the non-factorizing term is non-zero, and did not consider the specific case of QCD tree 
amplitudes. As mentioned in Section 2, it is the fact that the Feynman rules of QCD do 
not introduce any phases, so all of the elements of the colour vector Sp^\A4^) have the 
same phase, that leads to this result. In order to have a non-zero result, one must have 
a hard process with more than one colour flow, with non-trivial phase differences between 
them. Within the Standard Model, the only 2 — >■ 2 processes that satisfy these conditions 
are qq' — > qq 1 and qq — > qq in which an s-channel colour singlet propagator (W and Z 
respectively) introduces a non-trivial phase, which can then interfere with i-channel octet 
(gluon) or singlet (7 or Z) exchange. 



6 Conclusion and discussion 

We have explored the results presented in CdFR and rewritten them in a way that makes 
clear that the origin of the breakdown in collinear factorization is a mismatch in the non- 
Abelian colour matrix for Coulomb gluon exchange between the two incoming partons in 
the full matrix element and the sub-process matrix element it is factorized into, at one 
loop. At two- loop level, no fundamentally new non-factorizing effects enter, and the two- 
loop violation of strict factorization at the amplitude level is entirely determined by the 
one-loop one. This mechanism is identical to the one studied in Refs. [3, 4]. 
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At one-loop level the violation of factorization in the amplitude is anti-Hermitian and 
cannot contribute to physical cross sections. At two-loop level the sum of all contributions 
corresponding to two Coulomb gluon exchanges (in the amplitude, the conjugate amplitude 
and one in each) similarly cancels. However, a non-zero commutator between one Coulomb 
gluon exchange and one eikonal exchange gives a Hermitian term that may contribute to 
physical cross sections. We find that this contribution does vanish for pure QCD processes. 
Again, this mirrors the conclusions of Refs. [3, 4]. Ref. [4] in particular has a discussion of 
the cancellation of these terms in physical cross sections. 

In Refs. [3, 4], we concluded that the non-cancellation of Coulomb gluon effects does 
have a physical effect on non-inclusive QCD observables at the three-loop level. We showed 
that this comes about because only by exchanging two Coulomb gluons can we obtain a 
non-zero effect on physical cross sections and only via a non-zero commutator between 
these gluons and a further eikonal gluon can we avoid their complete cancellation as a pure 
phase effect. This then induced a mis-cancellation between the real and virtual initial- 
state collinear emissions (i.e. the breakdown of the 'plus-prescription'), which resulted in 
a "super-leading" logarithm. 

The all-order development of CdFR can be used at three loops to produce the results for 
the operator Iq , its non-factorizing part , and its effect on physical cross sections . 
Exactly the mechanism just described does indeed occur: the commutator between two 
non-factorizing one-loop Coulomb contributions and one factorizing eikonal contribution 
produces a non-factorizing term in the cross section that does not cancel. In order to 
display the results, we define one additional piece of notation: the Coulomb gluon parts of 



and i" are defined as A^ and A.p respectively (and therefore A^ 



A 



(i) 



A«). 

p ' 



i. c - ^ 

We obtain terms that are in one-to-one correspondence with those obtained for the first 
non-zero super-leading logarithm in Ref. [4]. In the 2-parton collinear limit, in which the 
eikonal part of Iq is a pure number, we obtain 



k (i) 

p ' 



A«7 (1) +7 (1) 

p ' 







A« 






p 






A« 






P 



p ' 



5p (0) . (6.1) 



The structure of the first line is identical to the contribution to super-leading logarithms 
from configurations in which the out-of-gap gluon has the highest k± (see Eq. (4.18) of 
Ref. [4]) and is driven by the mismatch between the double-commutator of one eikonal 
exchange and two Coulomb exchanges in the full and factorized matrix element. The 
structure and relative normalization of the second line are identical to the contribution 
to super-leading logarithms from configurations in which the out-of-gap gluon has the 
second-highest kj_ (see Eq. (4.20) 7 of Ref. [4]) and is driven by the mismatch between the 
commutator of one eikonal exchange and one Coulomb exchange in the full and factorized 
matrix element, with a second Coulomb exchange in the factorized matrix element giving 
an overall real Hermitian contribution. To obtain the result in the multi-parton collinear 
limit, one simply replaces J 1 by 1^ + I^q ■ 



7 This is Eq. (4.19) in the arxiv version of the paper. 
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We also note that the non- vanishing double-commutators in Eq. (6.1) are identical 
to those presented in Ref. [22] (Eq. (4.10)), where they triggered the failure of gluon 
Reggeization at NNLL. 

We should note that the results we have presented do not imply any mis-cancellation 
of the singularities of cross sections in hadron-hadron collisions, because at the longest dis- 
tance scales, which contribute to the infrared poles of the full hadron-hadron cross section, 
the incoming particles are the colourless hadrons and the colour factor associated with the 
factorization violation is zero. This was a crucial step in the proof of the factorization 
theorem in Ref. [1]: they were able to isolate the effect of the Coulomb/Glauber region of 
gluon momentum into factorization-violating terms that, when summed over all diagrams, 
gave zero. Furthermore, we fully expect that the factorization breaking effects we have 
been discussing will cancel in inclusive observables, in accord with the theorem derived 
by Collins, Soper and Sterman [1] and more recently by Aybat and Sterman [21]. The 
cancellation ought to occur upon accounting for real emissions with the same cut-diagram 
topology as the virtual corrections we consider here. However, these cancellations will 
not be complete in non-inclusive observables. In the case of super-leading logarithms, the 
factorization breaking is accompanied by an uncancelled soft-collinear double logarithm 
in an otherwise single logarithmic observable. The question of at which logarithmic order 
factorization breaking effects enter has been discussed by Banfi, Salam and Zanderighi [12]. 
In general, we expect to be able to factorize collinear logarithms only up to a scale \i be- 
low which real emissions are summed inclusively, and that the violations of factorization 
discussed here will contribute to the values of physical observables above that scale. 

Finally, we may remark that the factorization-breaking physics we have been exploring 
is not unrelated to the role of the underlying event, or of gap survival in diffractive physics. 
For an observable that is inclusive below some scale n, we may factorize the incoming par- 
tons at that scale. The evolution in scales between \jl and the high scale of a hard process 
will 'dress' these incoming partons with collinear accompaniments and interactions between 
these systems will violate the collinear factorization theorem, as we have discussed. These 
accompanying systems will interact with each other in an exactly analogous way to the 
remnants in hadron-hadron collisions. Thus, the factorization violations discussed in this 
paper will have effects that may be considered as perturbatively-calculable contributions to 
the underlying event. Furthermore, diffractive scattering is not power suppressed and can 
lead to important contributions in studies of the final state. However, it is well known that 
diffractive scattering in hadron-hadron collisions does not obey a factorization theorem: as 
a result one is led to the phenomenological idea of gap survival. In the high energy limit of 
scattering amplitudes, Coulomb gluon exchange is responsible for building the perturbative 
manifestation of the Pomeron and the eikonal exchange triggers the breakdown of factor- 
ization and leads to the filling of a would-be rapidity gap; in this case the factorization 
violating effects we have discussed can be thought of as the perturbative tail of the gap 
survival effect. 
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